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INTRODUCTION 

Light hadron spectroscopy Hes beyond the applicabihty regime of perturbative QCD. However, there is a rigorous and 
systematic expansion in the form of an effective field theory of QCD, known as Chiral Perturbation Theory (ChPT) [1], 
which provides a model independent description of the dynamics of the lightest mesons, namely, the Goldstone Bosons 
of the QCD spontaneous chiral symmetry breaking. Despite pure ChPT is limited to low energies and masses, here 
we review how, when combined with model independent dispersion relations, it leads to a successful description of 
meson dynamics, generating resonant states without a priori assumptions on their existence or nature. This " unitarized 
ChPT" is a useful tool to identify the spectroscopic nature of resonances through their dependence on the QCD number 
of colors Nc, but also to relate lattice results to physical resonances by studying their quark mass dependence. 

We will concentrate on the meson sector, where ChPT is most developed and converges somewhat better, and is 
built out of pion, kaon and eta fields only, as a low energy expansion of a Lagrangian respecting all QCD symmetries. 
Generically, it is organized in powers of 0{p'^/A^), where p stands either for derivatives, momenta or meson masses, 
and A = 4n fji, where f,t denotes the pion decay constant. ChPT is renormalized order by order by absorbing loop 
divergences in the renormalization of parameters of higher order counterterms, known as low energy constants (LECs) 
that carry no energy or mass dependence. Their values depend on the specific QCD dynamics, and have to be 
determined either from experiment or from QCD. The relevant remark for us is that, up to the desired order, the ChPT 
expansion provides a systematic and model independent description of how meson masses and amplitudes depend on 
QCD parameters like the light quark masses ih — (m„ + m^i)/! and m^, or the leading 1 /Nc behavior [2]. 

DISPERSION RELATIONS AND UNITARIZATION 

Elastic resonances appear as poles on the second Riemann sheet of the meson-meson scattering partial waves f/y of 
definite isospin / and angular momentum J. At physical values of s, elastic unitarity implies 

lmtij{s)^a{s)\tij{s)\^ ^ lm-^ = -a{s), f/y = — — I , with a{s)^2p/^fs, (1) 

ti]{s) Retjj-ia 

where s is the Mandelstam variable and p is the center of mass momentum. However, ChPT amplitudes, being an 
expansion tjj ~ fj^' + fjy^ H , with f"*^ = 0{p~''), can only satisfy Eq. (1) perturbatively 

lmtlj\s)^0, lmtlfis)^a{s)\t2is)\\ ... => lmtl'j\s)/tlf\s) = a{s), (2) 

and cannot generate poles. Despite the resonance region lies beyond the reach of standard ChPT, it can be reached 
combining ChPT with dispersion theory either for the amplitude [3] or for the inverse amplitude through the Inverse 
Amplitude Method (lAM) [4, 5, 6]. The first approach has been successfully used, combined with data on other 



channels and high energies, to, for instance, determine precisely the parameters of the /o(600) or k:(800) resonances. 
Unfortunately, this additional experimental input makes it difficult to relate these results to QCD parameters like A^^ 
or m. Hence we will concentrate on the one-channel JAM [4, 5], since it uses ChPT only up to a given order inside a 
dispersion relation - without additional input or further model dependent assumptions - providing an elastic unitary 
ampUtude with the correct ChPT expansion up to that order. Other unitarization techniques will be commented below. 

The one-loop ChPT Inverse Amplitude Method 

For a partial wave tu{s), we can write a dispersion relation (subtracted three times to suppress high energy 
contributions) 

Note we have explicitly written the integral over the right hand cut (or physical cut, extending from threshold, Sth to 
infinity) but we have abbreviated by LC the equivalent expression for the left cut (from to — oo). We could do similarly 
with other cuts, if present, as in the nK case. Similar expressions hold for and r^"*), but remembering that ^^^^ is a 
pure tree level amphtude and it does not have imaginary part nor cuts, they read: 

4^)=ao + ai., 4')=^o + M + ^2.V^ Im4^)(.0^.- (4) 

n Js,h -y {s' -s- le) " 

Note that from Eq.(l) the imaginary part of the inverse amplitude is exactly known in the elastic regime. We can then 
write a dispersion relation like that in (3) but now for the auxiliary function G = {t^*)^ /tij, i.e., 

where now PC stands for possible pole contributions in G coming from zeros in tij. It is now straightforward to expand 
the subtraction constants and use that Imfj^' = and Imfjy' = ajfjj'p, so that ImG = —\mt\^\ In addition, up to the 
given order, LC{G) ~ —LC{t^^), whereas PC is of higher order and can be neglected on a first stage. Then 



. . , r Imt^f is')ds' , f4-|, (2) (-41 



(5) 



We have thus arrived to the so-called lAM: 

tu^t^''/{t^'-t^), (6) 

that provides an elastic amplitude satisfying unitarity and has the correct low energy expansion of ChPT up to the order 
we have used. The PC contribution has been calculated explicitly [6] and shown to be, not just formally suppressed, 
but numerically neghgible except near the Adler zeros, away from the physical region. It is straightforward to extend 
the lAM to other elastic channels or higher orders [5]. Naively, by looking at (1), it seems that the lAM is derived by 
replacing ReJ^}' by its 0{p^) ChPT expansion. But, strictly speaking, (1) is only vahd in the real axis, whereas our 
derivation allows us to consider the amplitude in the complex plane, and, in particular, look for poles of the associated 
resonances. Let us remark that ChPT has been used always at low energies to evaluate parts of a dispersion relation, 
whose elastic unitarity cut is taken into account exactly. Thus, the lAM formula is reliable up to energies where 
inelasticities become important (even though ChPT does not converge at those energies) because ChPT is not being 
used there. Only when the energy is close to the Adler zero one should use a slightly modified version of the lAM [6]. 
When reexpanding, a few of the higher order terms are produced correctly by the unitarization but not the complete 
series — for a discussion of this issue for the scalar pion form factor see Ref.[7]. 

In Fig.l, we present some preliminary results [8] of an updated fit of the lAM nn and TzK scattering amplitudes 
to data, but simultaneously fitting the available lattice results on m;t,mK,fj[,fK and some scattering lengths. It is 
important to remark that the resulting LECs are in fairly good agreement with standard determinations: no fine tuning 
is required. As usual the fy (600), p (770), (^(SOO) and K* (892) appear as poles in the second Riemarm sheet of their 
corresponding partial wave. Actually, already ten years ago [5], with the elastic lAM we were able to generate poles 
for the p(770), ^r*(892) and the controversial a (or /o(600)), without any modeling of the integrals but just ChPT 
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FIGURE 1. Results of our updated lAM fit [8] (continuous line). We also show non-unitarized ChPT results with the LECs given 
in the second reference of [9] (dot-dashed line). . Two upper rows: lAM versus experimental data on and TiK scattering. Two 
lower rows: fit results compared to lattice calculations [10] of fji, /k, mn/fn and the Tt^Tt'^, K^K^, K^Tt^ scattering lengths. We 
fit up to rrin = 440MeV, but even beyond lattice results are well described (grey areas). Experimental references are detailed in [9]. 



approximations. The fact that resonances are not introduced by hand, but generated from first principles and data, is 
relevant because the existence and nature of scalar resonances is the subject of a long-lasting intense debate. The fact 
that the only input parameters are those of ChPT is very relevant because we then know how to relate our amplitudes 
to QCD parameters like Nc or the quark masses. 

Other unitarization techniques within the coupled channel formalism 

Naively one can arrive to (6) in a matrix form, ensuring coupled channel unitarity, just by expanding the real part of 
the inverse T matrix. Unfortunately, there is still no dispersive derivation including a left cut /or the coupled channel 



case. Being much more complicated, different approximations to ReT"^ have been used: 

• The fully renormalized one-loop ChPT calculation of Re ' provides the correct ChPT expansion in all channels, 
also with left cuts approximated to O(p^) [11,9]. Indeed, using LECs consistent with previous determinations within 
standard ChPT, it was possible [9] to describe below 1.2 GeV all the scattering channels of two body states made of 
pions, kaons or etas. Simultaneously, this approach [9] generates poles associated to the p(770) and ^*(892) vector 
mesons, together with the /o(980), ao(980), /o(600) and k (or/ro(800)) scalar resonances. 

• Originally [12], the coupled channel lAM was used neglecting the crossed loops and tadpoles. This approach is 
considerably simpler, and although it is true that the left cut is absent, its numerical influence was shown to be rather 
small, since the meson-meson data are nicely described with very reasonable chiral parameters and generates all the 
poles enumerated above. Let us remark that this approximation keeps the s-channel loops but also the tree level up to 
O(p^), and that this tree level encodes the effect of heavier resonances, like the rho. Thus, contrary to some common 
belief, this approach still incorporates, for instance, the low energy effects of t-channel rho exchange. 

• Finally, if only scalar meson-meson scattering is of interest, it is possible to use just one cutoff (or a dimensional 
regularization scale or a subtraction constant) that numerically mimics the combination of chiral parameters that appear 
in those scalar channels. This method - known as the '"chiral unitary approach'"- has become very popular, even 
beyond the meson-meson interaction realm, due to its great simplicity but remarkable success [13] and also because it 
is rather simple to relate to the Bethe-Salpeter formahsm [14] that provides additional physical insight on unitarization. 

With this method it has been shown [15] that, in the SU(3) limit, and assuming no quark mass dependence on the 
cutoff, all light scalar resonances degenerate into an octet and a singlet. 

Also with this method, but using a chiral Lagrangian for the pseudoscalar- vector interaction it has been possible to 
generate axial- vector mesons [16]. 

THE NATURE OF RESONANCES FROM THEIR LEADING 1 jNc BEHAVIOR 

The QCD 1 /A'c expansion [2], valid in the whole energy region, provides a rigorous definition of qq bound states: their 
masses and widths behave as 0(1) and 0{\/Nc), respectively. The QCD leading 1/A'c behavior of fn and the LECs is 
well known, and ChPT amplitudes have no cutoffs or subtraction constants where spurious Nc dependences could hide. 
Hence, by scaling with A^^ the ChPT parameters in the JAM, the mass and width Nc dependence of the resonances has 
been determined to one and two loops [17, 18]. These are defined from the pole position as ^Js^^ = M— iY. However, 
a priori, one should be careful not to take Nc too large, and in particular to avoid the N- — > °° limit, because it is a 
weakly interacting limit. As shown above, the lAM relies on the fact that the exact elastic RC contribution dominates 
the dispersion relation. Since the lAM describes data and the resonances within, say, 10 to 20% errors, this means that 
at A'c = 3 the other contributions are not approximated badly. But meson loops, responsible for the RC, scale as 3/Nc 
whereas the inaccuracies due to the approximations scale partly as 0{1). Thus, we can estimate that those 10 to 20% 
errors at A'c = 3 may become 100% errors at, say A'c ~ 30 or A'c ~ 15, respectively. Hence we have never shown results 
[17, 18] beyond A^^ = 30, and even beyond A^^ ^ 15 they should be interpreted with care. Of course, in special cases 
the 1AM could still work for very large A^^, as it is has been shown for the vector channel [19]. But that is not the case 
for the scalar channel, which, if used for too large Nc, leads to inconsistencies [19] for some values of the LECs. 

Thus, Fig.2 shows the behavior of the p, K* and a masses and widths found in [17]. The p and K* neatly follow 
the expected behavior for a qq state: M ~ 1, F ~ 1 /Nc. The bands cover the uncertainty in /i ~ 0.5 — 1 GeV where the 
LECs are scaled with A^^- Note also in Fig.2 (top-right) that, for that set of LECs, outside this ji range the p meson 
starts deviating from a a behavior. Something similar occurs to the K* (892). Consequently, we cannot apply the Nc 
scaling at an arbitrary jU value, if the well established p and K* qq nature is to be reproduced. 

In contrast, the a shows a different behavior from that of a pure qq: near Nc=3 both its mass and width grow with 
A'c, i.e. its pole moves away from the real axis. Of course, far from A'c = 3, and for some choices of LECs and n, the 
sigma pole might turn back to the real axis [18, 19, 20], as seen in Fig.2 (bottom-right). But, as coiimiented above, 
the lAM is less reliable for large Nc, and at most this behavior only suggests that there might be a subdominant qq 
component [18]. In addition, we have to ensure that the LECs used fit data and reproduce the vector qq behavior. 

Since loops are important in determining the scalar pole position, but are 1 /Nc suppressed compared to tree level 
terms with LECs, it is relevant to check the 0{p'^) results with an 0{p^) 1AM calculation. This was done within SU{2) 
ChPT in [18]. We defined a jj^^-like function to measure how close a resonance is from a qq Nc behavior. First, we 
used that jtf'^-like function at 0{p^) to show that it is not possible for the a to behave predominantly as a qq while 
describing simultaneously the data and the p qq behavior, thus confirming the robustness of the conclusions for Nc 
close to 3. Next, we obtained a 0{p^) data fit - where the p qq behavior was imposed - whose A'c behavior for the 
p and a mass and width is shown in Fig.3 (left and center). Note that both M„ and grow with A'c near A'c = 3, 
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FIGURE 2. Top: (Left and center) Nc behavior of the p and K* mass and width. (Right) Different p pole trajectories for different 
values of jx, note that for jx = \.2 GeV the p pole goes away the real sixis. Bottom: (Left and center) Nc behavior of the c mass 
and width. (Right) Different a pole trajectories for different |t values. 
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FIGURE 3. Left and center; Nc behavior of the p and <J pole at 0{p^) with the "p as qq fit". Right: Sigma behavior with Nc at 

0{p^) with the "p and a as qq fit". 



confirming the 0{p'^) result of a non qq dominant component. However, as Nc grows further, between Nc S and 
Nc ~ 15, where we still trust the lAM results. Ma becomes constant and Fa starts decreasing. This may hint to a 
subdominant qq component, arising as loop diagrams become suppressed when Nc grows. Finally, and despite this 
scenario is disfavored since the p starts deviating from its qq behavior, we checked how big this a qq component can 
be made. Thus we forced the a to behave as a qq using the above mentioned j^-like measure. We found that in the 
best case - Fig. 3. (right) - this subdominant qq component could become dominant around A'c > 6 — 8, at best, but 
always with an A'c ^ <» mass above roughly 1 GeV instead of its physical ~ 450 MeV value. 

Let us emphasize again [21] what can and what cannot be concluded from our results and clarify some frequent 
questions and doubts raised in this and other meetings, private discussions and the literature: 

•The dominant component of the a and k in meson-meson scattering does not behave as a qq. Why "dominant"? 
Because, most likely, scalars are a mixture of different states. If the qq was dominant, they would behave as the p or 
the K* in Fig. 2. But a smaller fraction ofqq cannot be excluded and is somewhat favored in our 0{p^) analysis [18]. 

•Two meson and some tetraquark states [22] have a consistent "qualitative" behavior, i.e., both disappear in the 
meson-meson scattering continuum as A'c increases. Our results are not able yet to establish the nature of that 
dominant component. To do so other tools might be necessary as, for instance, those outlined in [26, 27]. The most 
we could state is that the behavior of two-meson states or some tetraquarks might be qualitatively consistent. 

The Nc ^ '>° limit has been studied in [20, 19]. Apart from its mathematical interest, it could have some physical 
relevance if the data and the large Nc uncertainty on the choice of scale were more accurate. Nevertheless: 

• As commented above, a priori the lAM is not reliable in the Nc^°° limit, since it corresponds to a weakly interacting 
theory, where exact unitarity becomes less relevant in confront of other approximations made in the 1AM derivation. 
It has been shown [19] that it might work well in that limit in the vector channel of QCD but not in the scalar channel. 

• Another reason to limit ourselves to A'c not too far from 3 is that i n our calculations we have not included the 
Tj'(980), whose mass is related to the Ua{1) anomaly and scales as s/S/Nc. Nevertheless, if in our calculations we 
keep A'c < 30, its mass would be > 310MeV and thus pions are still the only relevant degrees of freedom for the 
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FIGURE 4. Left: Movement of the a (dashed hnes) and p (dotted hnes) poles for increasing m-n (direction indicated by the 
arrows) on the second sheet. The filled (open) boxes denote the pole positions for the CT (p) at pion masses m-ji = 1,2, and 
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respectively. For m-ji = 3m^'^^ three poles accumulate in the plot very near the nn threshold. Note that all poles are 
always far enough from the Adler zero (circles). Right: Comparison of our results for the Mp dependence on nijc with some recent 
lattice results from [30]. The grey band covers only the error coming from the LECs uncertainties. 



scalar channel in the a region. 
• Contrary to the leading 1 /Nc behavior in the vicinity ofNc — 3, the Nc — > °° limit does not give information on the 
"dominant component" of light scalars. The reason was commented above: In contrast to qq states, that become 
bound, two-meson and some tetraquark states dissolve in the continuum as Nc — > °°. Thus, even if we started with an 
infinitesimal qq component in a resonance, for a sufficiently large A^e it may become dominant, and beyond that A^e 
the associated pole would behave as a qq state. Also, since the mixings of different components could change with 
Nc, a too large A^^ could alter significantly the original mixings. 

Actually, this is what happens for the one-loop 1AM a resonance for A^^ °°, but it does not necessarily mean that 
the "correct interpretation [...] is that the cj pole is a conventional qq meson environed by heavy pion clouds" [20]. 
That the scalars are not conventional, is simply seen by comparing them in Figs. 1 and 2 with the "conventional" p and 
K* in those very same figures. A large two-meson component is consistent, but the A^^ °o of the one-loop unitarized 
ChPT pole in the scalar channel limit is not unique [20, 19] given the uncertainty in the chiral parameters. Moreover, 
for some LECs the scalar channel one-loop 1AM in the A^^- — > °o limit can lead to phenomenological inconsistencies 
[19], since poles can even move to negative squared mass values (weird), to infinity or to a positive mass square. That is 
one of the reasons why in the figures here and in [17, 18] we only plot up to N- = 30, but not 100, or a million. Hence, 
robust conclusions on the dominant light scalar component can be obtained not too far from real life, say A^^ < 15 or 
30, for a jj. choice between roughly 0.5 and 1 GeV, that simultaneously ensures the qq dependence for the p and K* 
mesons. Note, however, that under these same conditions the two-loop 1AM still finds, not only a dominant non-^<7 
component, but also a hint of a qq subdominant component, which is not conventional in the sense that it appears at a 
much higher mass than the physical a. This may support the existence of a second qq scalar octet above 1 GeV [23]. 

Finally, using not the 1AM, but the chiral unitary approach with a natural range for the cutoff A^^ dependence, it has 
also been suggested [24] that a large, in some cases dominant, non qq behavior could exist in axial vector mesons. 

QUARK MASS DEPENDENCE OF RESONANCES 

ChPT provides a rigorous expansion of meson masses in terms of quark masses (at leading order M^^^^^,^ ^ mg). Thus, 
by changing the meson masses in the amplitudes, we see how the poles generated with the 1AM depend on quark 
masses. In [25] we presented the SU(2) analysis for the p and a but here we also report on our recent developments 
[8] in the SU(3) formalism and the K'(800) and K*{892) strange resonances. 

The values of mjc considered should fall within the ChPT range of applicability and allow for some elastic nn and 
TtK regime below KK or Krj thresholds, respectively. Both criteria are satisfied if nij^ < 440 MeV, since SU{3) ChPT 
still works with such kaon masses, and because for m„ ~ 440 MeV, the kaon mass becomes ~ 600 MeV. Of course, 
we expect higher order corrections, which are not considered here, to become more relevant as m^^ is increased. Thus, 
our results become less reliable as nij^ increases due to the 0{p^) corrections which we have neglected 




FIGURE 5. ntj: dependence of resonance masses and widths in units of the physical values. In the two left panels the dark (light) 
band shows the results for the CT (p). The width of the bands reflects the uncertainties induced only from the uncertainties in the 
SU(2) LECs. Similarly, in the two right panels we show our recent developments [8] for the K"(800) and A^*(892) using central 
values of the SU(3) lAM fits. The dot-dashed (dotted) line shows the m,: dependence of the corresponding scalar (vector) width 
from the change of phase space only, assuming a constant coupling of the resonance to two mesons. 



Fig. 4 (left) shows the evolution of the a and p pole positions as nijc is increased. In order to see the pole movements 
relative to the two pion threshold, which is also increasing, all quantities are given in units of m„, so the threshold is 
fixed at y/s ~ 2. Both poles move closer to threshold and they approach the real axis. The p poles reach the real axis 
as the same time that they cross threshold. One of them jumps into the first sheet and stays below threshold in the real 
axis as a bound state, while its conjugate partner remains on the second sheet practically at the very same position as 
the one in the first. In contrast, the a poles go below threshold with a finite imaginary part before they meet in the 
real axis, still on the second sheet, becoming virtual states. As m„ is increased further, one of the poles moves toward 
threshold and jumps through the branch point to the first sheet staying in the real axis below threshold, very close to 
it as ruji keeps growing. The other a pole moves down in energies further from threshold and remains on the second 
sheet. These very asymmetric poles could be a signal of a prominent molecular component [26, 27], at least for large 
pion masses. Similar movements have been found within quark models [28] and a finite density analysis [29]. 

Fig. 4 (right) shows our results for the p mass dependence on m^j compared with some recent lattice results [30], and 
the PDG value for the p mass. Now the mass is defined as the point where the phase shift crosses 7z/2, except for those 
m„ values where the p becomes a bound state, where it is defined again from the pole position. Taking into account the 
incompatibilities within errors between different lattice collaborations, we find a qualitative good agreement with the 
lattice results. Also, we have to consider that the m„ dependence in our approach is correct only up to NLO in ChPT, 
and we expect higher order corrections to be important for large pion masses. The Mp dependence on m,: agrees nicely 
with the estimations for the two first coefficients of its chiral expansion [31]. 

In Fig. 5 (left) we compare the dependence of Mp and M(, (defined from the pole position .^i^^ = M — iT/2), 
normalized to their physical values. The bands cover the LECs uncertainties. We see that both masses grow with rriji, 
but Ma grows faster than Mp. Below m^n ~ 330 MeV we only show one line because the two conjugate a poles have 
the same mass. Above 330 MeV, these two poles lie on the real axis with two different masses. The heavier pole goes 
towards threshold and around m^j ~ 465 MeV moves into the first sheet, but that is beyond our appHcability limit. 

In the next panel of Fig. 5 we compare the niji dependence of Yp and normalized to their physical values: note 
that both widths become smaller. We compare this decrease with the expected phase space reduction as resonances 
approach the nn threshold. We find that Yp follows very well this expected behavior, which implies that the pnn 
coupling is almost m^i independent. In contrast, Ya deviates from the phase space reduction expectation. This suggests 
a strong m^i dependence of the a coupling to two pions, necessarily present for molecular states [27, 32]. 

Finally, in the last two panels of Fig. 5 we compare the mass and width dependence on m of the fc(800) versus the 
K*{%92), keeping fixed. Note that the same pattern of the a — p system is repeated. Belonging to the same octet, 
the K* (892) and p behave very similarly, and both their widths follow just phase space reduction. The a K behavior 
are only qualitatively similar, the latter being somewhat softer Among other effects, this might be due to a possible 
significant admixture of singlet state in the a. 

SUMMARY 



We have reviewed how the Inverse Amplitude Method (lAM) [6] is derived from the first principles of analyticity, 
unitarity, and Chiral Perturbation Theory (ChPT) at low energies. It is able to generate, as poles in the amplitudes, the 



light resonances appearing in meson-meson elastic scattering, without any a priori assumptions. Up to a given order in 
ChPT, it yields the correct dependences on the quark masses and the number of colors. 

The resonance leading 1/A^c behavior suggests that the dominant component of light scalars does not behave as a qq 
state as Nc increases not far from Nc = 3. When using the two loop lAM result in SU(2), below Nc ~ 15 or 30, there is 
a hint of a subdominant qq component, but arising at roughly twice the mass of the physical G. 

We have also predicted the evolution of the /o (600) , p (770) , (^(SOO) and K* (892) pole positions with increasing pion 
(quark) mass [25, 8] and have seen how they become bound states: softly in the vector case and with a non-analyticity 
in the scalar case. We have also shown that the vector-meson-meson coupling constant is almost m^i independent and 
we have found a qualitative agreement with some lattice results for the p mass evolution with nijt. These findings might 
be relevant for studies of the meson spectrum and form factors — see Ref . [33] — on the lattice. Work is in progress [8] 
to study also the strange quark mass dependence. 
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